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ABSTRACT 



We have combined 2MASS and POSS 11 data in a search for nearby ultracool (later than M6.5) 
dwarfs with Ks < 12. Spectroscopic follow-up observations identify 53 M7 to M9.5 dwarfs and seven 
L dwarfs. The observed space density is 0.0045 ± 0.0008 M8-M9.5 dwarfs per cubic parsec, without 
accounting for biases, consistent with a mass function that is smooth across the stellar/substellar limit. 
We show the observed frequency of Ha emission peaks at ~ 100% for M7 dwarfs and then decreases for 
cooler dwarfs. In absolute terms, however, as measured by the ratio of Ifa to bolometric luminosity, 
none of the ultracool M dwarfs can be considered very active compared to earlier M dwarfs, and we 
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show that the decrease that begins at spectral type M6 continues to the latest L dwarfs. We find that 
flaring is common among the coolest M dwarfs and estimate the frequency of flaxes at 7% or higher. 
We show that the kinematics of relatively active (EWhc > 6 A) ultracool M dwarfs are consistent 
with an ordinary old disk stellar population, while the kiucniatics of inactive ultracool M dwarfs are 
more typical of a 0.5 Gyr old population. The early L dwarfs in the sample have kinematics consistent 
with old ages, suggesting that the hydrogen burning limit is near spectral types L2-L4. We use the 
available data on M and L dwarfs to show that chromospheric activity drops with decreasing mass and 
temperature, and that at a given (M8 or later) spectral type, the younger fleld (brown) dwarfs are 
less active than many of the older, more massive fleld stellar dwarfs. Thus, contrary to the well-known 
stellar age-activity relationship, low activity in fleld ultracool dwarfs can be an indication of comparative 
youth and substellar mass. 

Subject headings: solar neighborhood — stars: activity — stars: kinematics — stars: 
low-mass, brown dwarfs — stars: luminosity function, mass function 

1. Introduction 

Catalogs of nearby stars (Gliese &; Jahreifi 1991; Kirkpatrick, Henry, &; Simons 1995; Reid, 
Hawley, &; Gizis 1995) and high proper motion stars (Luyten 1979) are grossly deficient in very 

low mass (VLM) dwarfs. With spectral types of M7 and later, these objects, sometimes called 
"ultracool M dwarfs,"^ are so optically faint that even nearby ones eluded searches based on the 
older (pre-1980s) sky surveys. These dwarfs have particular importance because they lie at or below 
the hydrogen burning limit - and have proven not only to be estimators of the numbers of dark 
brown dwarfs, but also present interesting astrophysical challenges in their own right. 

The new generation of sky surveys allows this deficiency to be addressed and large samples 
of nearby VLM dwarfs to be identified. The Two Micron All-Sky Survey^ (Skrutskie et al., in 
prep.; hereafter 2MASS) provides reliable photometry in the JHKg passbands, close to the peak of 
emission for these cool dwarfs. Furthermore, the Second Palomar Sky Survey (Reid et al. 1991), 
hereafter POSS II) provides Bj,Rf, and I^v photographic photometry in the northern hemisphere. 
In the southern hemisphere, the UK Schmidt and ESO sky survey plates provide Bj and Rp 
magnitudes. In sum, it is becoming possible to identify both the least luminous stars and young 
massive brown dwarfs by their optical and near-infrared colors alone over most of the sky. 

We present first results of a search using near-infrared and optical sky survey data aimed at 
completing the nearby star catalog for the ultracool M dwarfs. We discuss the sample selection and 
spectroscopic followup in Section 2. Although the sample discussed in this paper includes only a 
small fraction of the total population of nearby ultracool dwarfs, it represents a fourfold increase in 



^AU spectral types in this paper are on the Kirkpatrick, Henry, & McCarthy (1991) M dwarf and Kirkpatrick et 
al. (1999b) L dwarf systems. L dwarfs are cooler than "ultracool" M dwarfs. 

^2MASS data and documentation are available at http://www.ipac.caltech.edu/2mass 
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the number of such sources known. We discuss some preliminary results concerning the statistical 

properties of these sources in the latter sections of this paper. We discuss some stars of special 
interest in Section 3 and the 2MASS colors of ultracool M dwarfs in Section 4. The local space 
density of VLM dwarfs is discussed in Section 5, their activity and kinematics are discussed in 
Section 6, and finally our conclusions and future prospects are discussed in 7. 

2. Data 

2.1. Sample Selection 

Our results arc based on three observational samples. For our initial observing run, in July 
1998, wc used both photometric and proper motion criteria to define a sample of candidate VLM 
dwarfs. Based on the results from this run and further experience analyzing 2MASS data, we 
were able to improve our selection criteria for our December 1998 and subsequent observing runs. 
There are thus three samples with different properties, and when necessary we distinguish them as 
"Sample A" (July 1998), "Sample B" (December 1998), and "Sample C" (June 1999) respectively. 
Samples B and C have nearly identical selection criteria, and when combined are refered to as 
Sample BC. In all samples, objects within 20 degrees of the Galactic plane were excluded. Unless 
otherwise stated, all the analysis in Section 5 and 6 is based on Sample BC. 

Sample A was based on 2MASS data processed by July 1998. All these data were obtained 
at the Mt. Hopkins 2MASS telescope and a total of 363 square degrees were searched. All objects 
classified as extended by the 2MASS pipeline were eliminated (Jarrett et al. (2000)). The 2MASS 
data were correlated with the USNO's PMM scans of the POSS II plates, using a preliminary version 
of the software used by the 2MASS Rare Objects Core Project (this software will be described in 
more detail in a future publication by Monet et al.). This provided three additional colors: Bj, 
Hp, and Ijsf. Zero-point calibrations for the POSS II scans were not available, but we have since 
found that rough zero-points (good to ±0.5 magnitudes) are B = Bjinst + 2, Rc = Rpinst — 1) and 
Ic = iNinst — 1- These do not account for plate-to-plate variations or the significant color terms 
expected in Bj and Rp- 

All objects that met the following criteria were observed: 

1. 9.0 < < 13.0 

2. 0.95 < J-K^ < 1.30 

3. no Bj detection 

4. Significant (> 2a) proper motion 



Objects were also observed if they satisfied: 
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1. 9 < Ks < 12.0 

2. 0.95 < J-K,. < 1.30 

3. Bj detection 

4. Significant (> 2a) proper motion 

The proper motion criterion requires elaboration. The magnitude of the observed positional 
offset between the 2MASS and POSS II source is compared to the distribution of all 2MASS-P0SS 
II correlations. Only objects with significant proper motion were selected. Thus, while we selected 
sources that are moving with high confidence (roughly 2a), the actual cutoff in terms of arcseconds 
per year depends on the epoch difference between the F plate and 2MASS, which varies between 
and ~ 10 years. 

As shown in Section 2.2, these criteria led to the identification of a number of new nearby 
M dwarfs, but they are flawed in some respects. First, the POSS II and UKST Ilia J plates are 
sufficiently sensitive that many nearby (bright) VLM dwarfs arc detected. Moreover, the blue 
cutoff in J-Ks allowed mid-M dwarfs to enter the sample. We selected this blue cutoff because the 
prototypical M7 dwarf VB 8 has J-K= 0.95 (Leggett 1992). However, M dwarfs between spectral 
type MO and M6.5 all have J-K^ 0.9, and, with uncertainties of aj-x ~ 0.04 mag, significant 
numbers are scattered to J-K^ > 0.95. Consequently, our initial attempts to select ultracool M 
dwarfs produced a sample which is heavily contaminated by distant mid-M dwarfs. As it turned 
out, there were no 2MASS sources meeting our color and magnitude cuts that were not paired with 
a POSS II source in Sample A. 

Based on experience gained from this analysis, we revised our sample selection for the December 
1998 observing run. Because a larger area was available, we focused on brighter VLM dwarfs which 
should lie within ~ 20 parsecs. The Sample B criteria are: 



1. 


Ks < 12.0 


2. 


J-K^ > 1.00 


3. 


Rf -Ks> 3.5 or I;v - > 2.0 


4. 


S < +30deg 


5. 


a < 13^^00"" or a > 20^00"* 


6. 


J-H < 1 H-K +0.25 



No selection based on proper motion was applied, and therefore Sample B is kinematically 
unbiased. A total of 2977 sq. degrees were searched. The J-K^ cutoff excludes early and mid-M 
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dwarfs, but also some of the bluer M7 dwarfs. Nearby bright L dwarfs, however, are included in 
these samples, since we impose no red cutoff. Both samples arc magnitude selected, and therefore 
arc biased towards ovcrluminous stars and unresolved near-equal luminosity binaries, although 
binaries with separations of a few arcseconds may be excluded by the extended source provision. 
The position selection is due to the requirement that the objects be observable from Las Campanas 
in December. Ultracool M dwarfs have R-K^ > 5.5 and I-K> 3.4 (Leggett 1992), but we used a 
more liberal selection to allow for uncertainties in the calibration of the photographic magnitudes. 
In Section 4, we show that R-K^ > 4.9 includes all M8 and later dwarfs. The J-H,H-K cut excludes 
M giants. 

Sample C was selected for our June 1999 Kitt Peak observing run. The selection was identical to 
the Sample B selection, except that a different area was covered and only objects with R — Kg > 4.9 
were selected. Processed data which lay within the following limits were selected: 

1. a > ll'^OO"* and 5 > +6° 

2. 16^^20"^ <a< 23^^35™ and -36 < S < +6° 

2.2. Spectroscopy and Data Analysis 

Sample A was observed on UT dates July 30 - August 1 1998 using the Double Spectrograph 
and the Hale 200-in. telescope during an observing run that was primarily devoted to our ongoing 
spectroscopic survey of Luyten high proper motion stars (Gizis Sz Reid 1997). The wavelength 
coverage included 6290 to 8800 A at a resolution of 3 A. 

Sample B was observed on December 2 - 7 1998 using the Modular Spectrograph on the Las 
Campanas 100-in telescope. The "Tek 5" chip, a 2048-square CCD with 24/x pixels, was used with 

a 600 1/mm grating blazed at 7500 A. The useful wavelength range of the spectra was 6100 - 9400 
A at a resolution of 6 A. A few targets (including all three objects with J-K^ > 1.3, which had been 
previously identified as L dwarf candidates) were observed during Keck observing runs (see Kirk- 
patrick et al. 1999b, hereafter K99, and Kirkpatrick et al. 2000, hereafter KOO) using LRIS (Oke et 
al. 1995). The resolution was 9 A with wavelength coverage from 6300 to 10100 A. Observations of 
the flare dwarf 2MASSI J0149089+295613 have already been described in Liebert et al. (1999). The 
Ha activity levels adopted here are the average quiescent values. The known object BRI 1222-1222 
was not observed, and we rely on the spectroscopic observations reported by Kirkpatrick, Henry, 
& Simons (1995) and Tinney & Reid (1998, hereafter TR). 2MASSW J0354013-h231633 has been 
previously published as 2MASP J0354012+231635 (Kirkpatrick, Beichman, & Skrutskie 1997) but 
the spectral observations presented here are independent. 

Sample C was observed June 22 - 23 1999 using the R.C. spectrograph and the Kitt Peak 4m 
telescope. The wavelength coverage was 6140 A to 9200 A using the 2048 CCD, but the extreme 
ends of the spectra were out of focus. A few objects were observed at Palomar 200-in. in May 1999. 
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The four new L dwarfs identified with the Kitt Peak data were reobserved at higher signal-to-noise 
using Keck in July 1998. The spectral measurements used for classification on the K99 system are 
given in Table 2. 

All spectra were extracted and flux calibrated using IRAF. M dwarf spectral types were mea- 
sured by overplotting dwarfs of known spectral type, and should be good to ±0.5 subclasses. All L 
dwarfs have Keck observations and were classified as in K99 and KOO. A few M dwarfs in Table 1 
have classifications that differ by 0.5 subclasses from previously published values - we leave our 
values unalterred as representative of our uncertainties. Ha fluxes were measured assuming the 
data were photometric (there were no clouds for our observations in December 1998). We assume 
that BCk = 3.2 as derived by Tinney, Mould, & Reid (1993) for ultracool M dwarfs. We believe 
the Ha fluxes should be viewed with caution since slit losses and the high airmass of many of the 
observations will increase the uncertainties (the spectrograph was not adjusted to the parallactic 
angle). However, our derived Ha to bolometric luminosity ratios are consistent with those of TR, 
and in any case the Ha emission strength is variable in these dwarfs. In Table 3, we list our 
measured and derived parameters for the ultracool dwarfs in Sample BC. 

Proper motions were estimated by measuring positions off the DSS (POSS I/UKST J) and 
XDSS (POSS II/UKST R) images of the photographic sky surveys. When a second epoch photo- 
graphic sky survey was not available, we used the 2MASS images instead. All motions reported are 
relative to other stars in the field, but the correction to absolute motions is negligible compared to 
other sources of error in the kinematics. The proper motion reported for LHS2397a is from Luyten 
(1979) and BRI1222-1222 is from Tinney (1996). Note that all of our targets are visible on the DSS 
images, but most were previously unrecognized. The targets which had no POSS H pairings in our 
initial processing all proved to have high proper motions. They are visible on the XDSS but lie 
outside the 8 arcsecond search radius employed in cross-referencing against the photographic data. 
Since the 2MASS positions are highly accurate, and both the 2MASS images and DSS images are 
(or will shortly be) easily accessible electronically, we are not presenting finding charts. 

Using the available parallaxes for late-type dwarfs (Monet et al. 1992; Tinney et al. 1995; 
Tinney 1996; Kirkpatrick et al. 1999b), we find (Figure 1) the linear fit Mk = 7.593 + 2.25 x 
J-Ks This fit is only valid for M7 and later dwarfs over the color range 1.0 < J-K^ < 1.6, and 
should be modified as more L dwarf parallaxes are measured at USNO. The observed scatter is 
a = 0.36 magnitudes. The distances and tangential velocities derived using this estimate are listed 
in Table 3. We caution, however, that the distances derived for many of the M7 and M7.5 dwarfs 
may be underestimated in this paper. As can be seen in Figure 1, the main sequence bends sharply 
at spectral type M7. Our selection of only targets with J-K^ > 1.0 tends to select M6's and M7's 
with overestimated colors, while a spectral classification error of only 0.5 subclasses from M6.5 to 
M7.0 leads to a large error in Mk- 
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3. Stars of Special Interest 

A few of the targets deserve special comment. Our search has identified one very nearby star 
and a number of very high proper motion dwarfs. The M8.0 dwarf 2MASSW J0027559+221932 has 
a photometric parallax that places it within ten par sees; given the uncertainties, it may lie within 
the eight parsec sample. Wc note also that a number of dwarfs in Table 3 have motions greater 
than 1 arcsccond per year, but do not appear in the LHS Catalog even though they are visible on 
the POSS plates. 

Seven L dwarfs are part of our Sample BC. The L5 dwarf 2MASSW J1507476-162738, also 
in Reid et al. (2000), was selected for this project but docs not lie in the Sample BC area. Addi- 
tional observations and discussion of 2MASSI J0746425+200032, 2MASSW J0036159+182110 and 
2MASSW J1439283+192915 are given in Reid et al. (2000). 2MASSW J1439283+192915 is in the 
original K99 paper while 2MASSW J0036159+182110 was observed at Keck for the KOO paper. 
2MASSW J1300425+191235 is particularly surprising: it has the J-K^ color of an M8 dwarf but 
has an LI spectrum. It will be discussed further in an additional paper, but we note here that 
the estimated distance and tangential velocity is based on the J-K^ color and should be viewed 
with great caution. Nevertheless, it apparently has a high velocity, and is likely to be old. While 
we report photometric distance estimates only, most of these dwarfs are on the USNO parallax 
program and accurate distances should be forthcoming. 

Since our selection is based upon photometry only, we are sensitive to wide binary pairs 
where the 2MASS observations of the secondaries are unaffected by the primaries. Two M dwarf 
secondaries that do not meet the spatial restrictions and were specially observed have been reported 
in Gizis et al. (1999b). One of the Sample C ultracool M dwarfs also appears to be a secondary. 
The M8.0 dwarf 2MASSW J2331016-040618 is 447 arcseconds west and 65 arcscconds south of the 
F8 dwarf HD 221356. Our photometric parallax of 26.3 parsecs for the M dwarf is consistent with 
Hipparcos trigonometric parallax of 26.24 parcsecs (Ferryman et al. 1997), as are the observed 
proper motions. This is apparently a wide binary system with separation of 0.057 parsecs. The F8 
primary may provide a useful age and composition constraint on the M dwarf. 

Two of our sources have been previously identified as candidate Hyades members. LP 475-855 
has been discussed as an Hyades candidate by Eggen (1993) although it was rejected as too bright 
by Leggett, Harris, &: Dahn (1994). Our photometry supports the latter conclusion, although 
conceivably it could be a foreground (escaping?) Hyades member if it is an unresolved equal-mass 
binary. Our initial observation of this star found it in a flare state, with Ha equivalent width of 
40 A. A 25 December 1998 Keck spectrum found Ha of only 7 A, which is the value we report in 
Table 1. LP 415-20 (Bryja 262) has been extensively discussed as a Hyades member and has been 
classified as an M6.5 dwarf (Bryja, Humphreys, & Jones 1994). The difference in spectral types is 
within our uncertainties, and we note that VO is visible in Bryja's plot of the spectrum. The poor 
distance estimate for this object is consistent with our belief that our M7 distances based on J-K^ 
colors are unreliable and should be viewed with caution. 
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4. 



Colors 



In Figure 2, we plot the R-Ks,J-Ks diagram for Sample B. We concluded that adjusting our 
color criterion to Rf-K^ > 4.9 would increase our selection efficiency without losing MS and later 
dwarfs, and we adopted this selection criterion for Sample C. This observation is consistent with 
our estimate that our simple Rp zero-point is good to ±0.5 magnitudes. 

Our observations show a good correlation between the far-red spectral type and the 2MASS 
near-IR colors. In Figure 3, we plot the observed J-K^ color distribution as a function of spectral 

type. As expected, the M7.0 and M7.5 distributions are truncated by our requirement that J- 
Ks > 1.0. The histograms suggest that of order one M8-M8.5 dwarf may be expected to be missed 
due to this requirement, and it appears that it is very unlikely than any (normal) M9 dwarfs are 
missed. 2MASS photometry for these bright sources is expected to be good to 0.03 magnitudes. 

In Figure 4, we plot the 2MASS near-infrared color-color diagram for our sample. It is evident 
that the M/L dwarf sequence lies well below our imposed J-H,H-K cut, and therefore we are not 
missing ultracool M dwarfs due to this criteria. We fit the relation J — Kg = (0.146 it 0.117) -|- 
(1.238 dr 0.263) X H — Kg assuming the errors in each color are 0.042 and including the M7.0 to 
M9.5 dwarfs. This relation may be convenient as a representative sequence in the 2MASS color 
system. It is interesting to note that the dwarfs around (H-K,J-H) = (0.45, 0.62) are nearly all 
classified as M7-M7.5 dwarfs, while the dwarfs at (II-K,J-H) = (0.42,0.7) but with similar J-K^ 
colors are nearly all classified as M8-M8.5 dwarfs. This may reflect some relation between the red- 
optical region (dominated by TiO and VO, and influenced by dust) and the IR colors (dominated 
by H2O and H2, and also influenced by dust) - or it is due to some subtle bias in our classifications 
or photometry (for example, perhaps we tend to select M7 dwarfs whose H-K^ color has been 
overestimated). Note that one of the outliers below the normal J-H,H-Ks relation is the peculiar L 
dwarf 2MASSW J1300425+191235. 



Our Sample BC is the first large sample of bright, photometrically selected ultracool M dwarfs. 
Using our data and derived distances, we can estimate the luminosity function using Schmidt's 
(1968) V/Vmax technique. The space density is 



5. Luminosity Function 




V^n.a. = f (lO.O(^«--^-+^-°)/^-°)' 



In our case, Ku^ = 12.0 and U = 6040 sq. degrees. The corresponding variance is 




max 



-9- 



The space densities are given in Table 4. We derive a space density of 0.0045 ±0.0008 ultracool 
M dwarfs per cubic parsec. According to our adopted color-magnitude relation, this is for dwarfs 
in the range 9.8 < Mk < 10.8. Therefore, the corresponding luminosity function bin is <I>(Mx = 
10.3) = 0.0048 ± 0.0009 dwarfs per cubic parsec per K magnitude. Since Tinney, Mould, & Reid 
(1993) have shown that BCk = 3.2 for these late dwarfs, this may be represented in bolometric 
magnitudes as ^{Mjjoi = 13.5) = 0.0048 ±0.0009 dwarfs per cubic parsec per bolometric magnitude. 
We note, however, that this value excludes the M7 dwarfs which will also contribute near M;^- = 9.8, 
so our value is a lower limit. The space density for the early L dwarfs is half that of the ultracool 
M dwarfs, although we caution that the distance estimate for 2M1300 may be incorrect. 

Schmidt's statistic measures the uniformity of the density distribution of a sample, effectively 
providing an estimate of sample completeness. For a uniform sample, < V/Vmax >= 0.5 with an 
uncertainty of J— where N is the number of stars observed. Both our L dwarf and M8.0-M8.5 



sample lie within 1 cr of this value, suggesting that we are complete. The value for M9.0 to M9.5 
dwarfs is more problematic, indicating that we have either excluded a few nearby, very bright M9 
dwarfs, or that there happen to be no such very nearby dwarfs in our survey volume. 

Our space density for the ultracool M dwarfs is consistent with Tinney 's (1993) value of 
$(M(,oi = 13.5) = 0.0076 ± 0.0031 dwarfs per cubic parsec per bolometric magnitude, which was 
based on selection with Rp and In photographic magnitudes but K followup of all VLM dwarfs 
to improve photometric parallaxes. Only 6 dwarfs contributed to this bin, accounting for Tinney's 
larger uncertainty relative to our sample. Delfosse et al. (1999) have analyzed the DENIS Mini- 
survey and found 19 MS and later dwarfs, including 3 L dwarfs. They do not estimate the M dwarf 
space density, but use the three L dwarfs to estimate $(M{,o; = 15.3) > 0.011 ± 0.006 dwarfs per 
cubic parsec per bolometric magnitude. We note that their estimated Mk for the ultracool M 
dwarfs are inconsistent with our adopted values (and Figure 1) since they consider their M8-M9 
dwarfs to have Mk > 11. 

Malmquist bias will affect our sample. Stobie, Ishida &; Peacock (1989) have shown that the 
luminosity function will be overestimated by: 



A model luminosity function can be used to derive the first and second derivatives. The scatter 
of parallax stars about the linear fit adopted here is cr = 0.36. Adopting this value for a, and 
making the assumption that the luminosity function is flat, we find that ^ = 0.21: i.e., the 
values we derive overestimate the true space densities by ~ 20%. Since this is only a preliminary 
sample, we defer further analysis of the Malquist bias, as well as the effect of unresolved binaries, 
until additional data are available. In the long term, trigonometric parallaxes and searches for 
companions will allow this issue to be addressed directly. 



— = 0.5(7^ (0.6 In 10)^ - 1.2 In 10— ± — 



Our derived space density can best be compared to the luminosity function of nearby stars. 
Figure 5 plots the Reid & Gizis (1997) {5 > —30°) luminosity function for stars within eight parsecs 
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with our M8.0-L4.5 data point added to it. Our results suggest that the dropoff seen in for the 
faintest {Mk > 10) dwarfs in the eight parsec sample is in part due to incompleteness. ^ Applying 
standard mass-luminosity relations to ^(Mk) derived from the 8-parsec sample implies a turnover 
in the mass function close to the hydrogen-burning limit. There is, however, no reason to expect 
the star formation process to be cognizant of the mass limit for hydrogen burning. Reid et al. 
(1999a) have modelled a sample of twenty 2MASS and and three DENIS L dwarfs and conclude 
that the substellar mass function is consistent with an extension of the power-law matched to data 
for stars with masses between 0.1 and IMq. The higher space densities measured for ultracool 
dwarfs in this paper suggest a greater degree of continuity across the stellar/substellar boundary. 
Continuation of the present survey should identify the missing dwarfs within eight parsecs. 

If compared to the classical "photometric" luminosity functions (Stobie, Ishida & Peacock 1989; 
Tinney 1993) which have a peak at Mf,oi = 10 and a dropoff to M^oi = 12, then our data would imply 
a rise in the luminosity function at the stellar /substellar boundary. However, this peaJi and dropoff 
are an artifact of the data analysis due to the incorrect assumption of a linear color-magnitude 
relation (Reid &; Gizis 1997) and/or other systematic errors such as unresolved binaries (Kroupa, 
Tout & Gilmore 1993). We believe that the nearby-star sample is a better comparison for our 
sample, and we emphasize again that the M dwarfs here identified require follow-up trigonometric 
parallax determinations and high-resolution imaging and radial velocity searches for companions 
to produce a definitive luminosity function. 

6. Activity and Kinematics 
6.1. Review 

The BC sample was not selected on the basis of proper motions, and therefore is (relatively) 
unbiased in terms of kinematics.^ It is useful to review the properties of nearby disk stars and 
the already-known properties of ultracool M dwarfs before discussing our kinematic and activity 
measurements. 

Stars are born with low space velocity dispersions and and high chromospheric activity levels. 
Over time, the space velocity of the stars increases as they interact with the Galactic disk. Using 
BAFGK dwarfs with known ages, Wielen (1974) showed that the total space velocity increases from 
o'tot = 19 km/s at mean ages of 0.4 to 0.9 x 10^ yr to 34 km/s at 2 x 10^ yr to 48 km/s at 5 x 10^ 
yrs. The high chromospheric activity levels of young stars with convective envelopes is attributed 



''Note that, as seen in Reid & Gizis's Fig. 2, the oft-used 5.2 parsec sample shows the same feature, albeit with 
less significance due to the very small volume. 

®The existence of age-luminosity, age-metallicity, metallicity-luminosity, age-activity, and age-kinematics corre- 
lations implies that there may bo kinematic and/or activity bias due to our Malquist-typc luminosity bias. If the 
distances are underestimated due to bias, then the estimated tangential velocities will also be biased. 
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to a dynamo which is driven by rotation. As the star ages, angular momentum loss through the the 

stellar wind spins down the star, causing the chromospheric activity to in turn decrease. Wielen 
(1974) did indeed find that Ca II emission line strength in M dwarfs was related to kinematics, 
with older stars showing less activity and higher space velocities. As progressively less massive 
(later spectral type) stars are considered, the observed frequency of high Ha activity increases. 
This is not due only to the fact that Ha emission is more detectable against the cool photosphere - 
Hawley, Gizis, & Rcid (1996, hereafter HGR) showed that the percentage of highly active M dwarfs 
increases with cooler spectral types even when Ha activity is compared to the star's bolometric 
luminosity. The increased lifetime of activity is confirmed by observations of open clusters Hawley, 
Reid, & Tourtellot (1999). The connection to rotation in early M dwarfs is confirmed by Delfosse 
et al. (1998), who have found that the incidence of rapid rotators is higher among cooler M dwarfs, 
and that the rapid rotators are active. There is some evidence that the rotation-activity relation is 
breaking down (Hawley, Reid, &; Gizis 1999). In summary, mid-M dwarfs maintain H a emission 
for billions of years as they slowly spin down. Kinematics are a good age indicator, but only in a 
statistical sense - individual stars can not be accurately dated by the velocities. 

For the ultracool M dwarfs, there is considerable evidence that the standard stellar age-activity 
and rotation-activity relations no longer apply. Basri & Marcy (1995) found that the M9.5 dwarf 
BRI 0021-0214 had very rapid rotation (vsini = 40 km/s) but little Ha activity. TR have shown 
that the lithium M9 brown dwarf LP 944-20 (Tinney 1998) is also a member of this class of 
"inactive, rapid rotators" as are the two of the DENIS L dwarfs (Martin et al. 1997). TR argue 
that observations of these field objects as well as open clusters indicate that the violation of the age- 
activity connection is primarily correlated with mass (the physical mechanisms remain unknown). 
Basri (1999) reports that rapid rotation is more common among objects of lower luminosity, and 
proposes that the Ha activity is powered by a turbulent dynamo that is quenched at high rotation 
rates. There is some evidence that even among the "inactive, rapid rotators" log ^xft^ is related 
to age, as it decreases from ~ —4.6 for the Pleiades to ~ —5.2 for ~ 0.5 — 1.0 x 10^ yr brown dwarfs. 
This trend may be true for younger ages, since the very low mass (~ 0.01 — O.O6M0), very young 
(< 10 Myr) M8.5 brown dwarf p Oph 162349.8-242601 has EWna > 50A (Luhman, Liebert & 
Rieke 1997; Martin, Basri, & Zapatero Osorio 1999), while the young 1 Myr), possible brown 
dwarfs (M « O.O7M0) V410 Tau X3 (M8.5) and X6 (M6) have EWna « 15 A (Martin, Basri, 
&; Zapatero Osorio 1999). Both imply higher activity levels than their older Pleiades and field 
counterparts, though it should be noted that they are also lower mass. In contrast to the "inactive, 
rapid rotators'," some ultracool M dwarfs do show Ha emission, but they have lower rotation rates 
(< 20 km/s, TR). However, even a rotation rate of ~ 5 km/s is adequate to maintain Ha emission 
in mid-M dwarfs (Bopp & Fekel 1977), so the limits on the "low" rotation rates of these ultracool 
M dwarfs are not surprising by comparison. A so-far unique field object is the M9.5 dwarf PC 
0025-1-0447 (Schneider et al. 1991; Martin, Basri, & Zapatero Osorio 1999), whose quiescent Ha 
emission (300 A) is comparable to highly active mid-M dwarfs in terms of The nature of 

this object and its emission is uncertain — Martin, Basri, & Zapatero Osorio (1999) argue that 
this object is a very young brown dwarf, suggesting that the USNO parallax indicating ordinary 
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ultracool M dwarf luminosity is incorrect. Some of the known M8-M9 dwarfs are definitely brown 

dwarfs. The "inactive, rapid rotator" M9 dwarf LP 944-20 has lithium absorption and a luminosity 
that indicates it has a mass between 0.056 and 0.064 Mq (Tinney 1998). The Pleiades brown 
dwarfs Teide 1 and Calar 3 have spectral types of MS, lithium absorption, and anomalous VO and 
Na features due to low-surface gravity (Martin, Rebolo, & Zapatero-Osorio 1996). 

6.2. The Properties of M and L dwarfs 

The BC sample of ultracool M dwarfs provides the first opportunity for a thorough investigation 
of the distribution of activity in these VLM dwarfs. In Figure 6, we compare the percentage of 
ultracool M dwarfs observed in emission to the HGR statistics for nearby K7 to M6.5 dwarfs. Note 
that that 80 to ~ 300 dwarfs contribute to each of the HGR bins up to spectral type M4.5. Only 
a few objects contribute to the M6.0 and M6.5 bins, which also are probably kinematically biased 
against young, active stars due to incompleteness in the pCNS3.^ We extend our M dwarf sample 
to even cooler dwarfs by using the K99 and KOO data, who report the strength of Ha emission in 
their L dwarfs. Their sample is photometrically-selected and kinematically unbiased. The sample 
shows a steady decline in Ha emission frequency from 60% (80% if a marginal detection is included) 
for type LO down to only 8% (25% if two marginal detections are included) for L4 dwarfs. None 
of the twenty dwarfs with spectral type L5.0 or later show definite emission (two have marginal 
detections). For the earlier dwarfs, emission of 1 A equivalent width would have been detected in 
almost all the objects; however, the upper limit on the equivalent widths for the latest L dwarfs 
was typically somwhat larger. The fact that there is so little photospheric continuum for the latest 
L dwarfs around the Ha feature should compensate for the lower sensitivity in terms in equivalent 
widths. 

The data indicate that the frequency of emission increases with later spectral type (cooler 
temperatures), until at spectral type M7 all of our targets show detectable emission. Indeed, we 
are not aware of any inactive M7 dwarfs (HGR; Gizis & Reid 1997). This indicates that the dwarfs 
can maintain detectable levels of activity for the lifetime of the Galactic disk. Later than M7, the 
Ha emission frequency begins to decrease, with our sample of ultracool M dwarfs merging cleanly 
with the L dwarfs. This coincides with the breakdown of the rotation-activity relation already 
noted for M9 and L dwarfs and reflects the apparent relative inability of the ultracool M dwarfs 
to heat the chromosphere discussed by TR. The percentage of emission at spectral type M6 is 
particularly uncertain, as the HGR sample may be biased toward higher velocity, hence older, stars 
at such low luminosities. Sixteen of our nineteen M6-M6.5 dwarfs show emission, but they have 
been effectively selected on the basis of unusually red J — Kg colors, and may be biased in some 
way. In any case, there is little doubt that some high velocity, presumably very old, M6 dwarfs are 



HGR's statistics for M7 and later dwarfs are sparse, but further observations have revealed that they are incorrect. 
Their Table 5 should show that 2 of 2 M7 dwarfs, 2 of 2 M8 dwarfs, and 2 of 3 M9.0-M9.5 dwarfs show emission. 
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no longer active. 

Since tlie Ha line is seen against an increasingly faint photosphere for these ultracool dwarfs, 
the log (^xfy^ ratio is more indicative of the true level of activity. We plot our M dwarfs, the HGR 
early M dwarfs, and the K99 L dwarf data in Figure 7. For the K99 data, we have measured the 
continuum level off the observed spectra and assumed BCk = 3.33 as measured by Tinney et al. 
(1993) for GD 165B to convert the K99 equivalent widths to flux ratios. Note that the increase 
in maximum observed activity levels from K7 to the peak at M3-M5 reflects the increased lifetime 
of emission for the lower mass stars. Earlier M dwarfs with activity levels near -3 are known in 
young clusters, but do not have long enough lifetimes to appear in the local sample. The lower 
envelope of data points is set by the fact that the minimum observable Ha emission of lA equivalent 
width corresponds to a lower luminosity fraction in cooler dwarfs. The addition of our data to the 
HGR data clearly indicates that beyond M7 the level of activity is indeed declining. This decline 
continues for lower (L dwarf) temperatures. The decline is quite steep - in only three subclasses 
(M8 to LI) the activity drops by one full dex. Martin et al. (1999) note a "slight trend toward 
decreasing emission in the L dwarfs" - our conclusions differ due to our larger sample and their use 
of equivalent widths only. 

While the quiescent Ha chromospheric activity is declining, our data suggest that flare activity 
is common in the ultracool M dwarfs, and may be a significant contributor to the activity energy 
budget. We summarize evidence for variability in Table 5. Since these events represent a strong 
enhancement of the Ha line strength, we suggest that they may be flares. At least a few dwarfs 
apparently maintain strong quiescent emission - our Ha line strength of 29 A for LP412-31 is 
identical with the value observed by Martm, Rebolo, & Zapatero-Osorio (1996). Other ultracool 
M dwarfs not in our sample have been seen to flare - Reid et al. (1999b) recently observed a flare 
on the "inactive" dwarf BRI 0021-0214, while Martin, Rebolo, & Magazzu (1994) observed the 
H-alpha EW of LHS 2065 on two consecutive nights as 7.5 A then 20.3 A. Martin et al. (1999) 
have observed flares in a number of ultracool M dwarfs. Flaring activity has thus been observed in 
ultracool M dwarfs with both very low and high levels of quiescent emission. 

Assuming that these strong variations are due to flares, we estimate the flare rate from our own 
statistics. At least four of the fifty-three ultracool M targets were flaring the first time we observed 
them for this program — implying that ultracool M dwarfs spend > 7% of the time in a flare state. 
This is consistent with the observed flare rates of the "inactive" M9.5 dwarf BRI 0021-0214 (Reid 
et al. 1999b) and the monitoring of 2M0149 (Liebert et al. 1999). This flare rate is a lower limit, 
since some of our other targets may also be flaring, but lacking additional spectra we cannot tell 
whether they are merely active as for LP 412-31, and since we have no way of identifying weaker 
flares. The Ha equivalent widths appear to be enhanced by a factor of ~ 10 in the observed flares, 
implying that perhaps half of the Ha luminosity is emitted during flares. 

We now consider the relationship between kinematics and activity for the M8 - M9.5 dwarfs. 
In Figure 8, we plot the observed relation between the tangential velocity and Ha emission. There 
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is a striking relationship between activity and velocity, in a sense opposite to that observed in 
more massive M dwarfs. All the dwarfs with strong Ha emission have large velocities {vfan > 20 
km/s). There is also a striking population of low- velocity, low-activity ultracool M dwarfs. With 
the exception of the high-velocity, inactive M9 dwarf 2MASSW J0109216-h294925, the least active 
stars appear to be drawn from a lower velocity population. It is difficult to fairly characterize 
the tangential velocity dispersion (cr^^^ = Y''of^i~+~^^) of these populations, but the inactive, 
low- velocity population in the lower left of Figure 8 may be characterized by atan ~ 15 km/s. 
While the low-velocity, low-activity population seems to to have EWho ^ 7A and vtan ^ 25 
km/s, we can calculate a velocity dispersion only for a purely activity selected sample. As an 
illustration, the dwarfs with EWHa < SA (excluding 2MASSW J0109216-F294925) have atan = 13 
km/s; in contrast, those with more emission have atan = 38 km/s. Using the approximation that 
atot = \J\crtant the implied total space dispersions of the two populations are 16 km/s and 47 
km/s. Comparing to Wielen (1974), the active M dwarfs are apparently drawn from a ~ 5 x 10^ yr 
population, but the inactive M dwarfs axe consistent with a ^ 0.5 Gyr population. This estimate 
is crude at best, but it seems clear that the overall ultracool M dwarf population is drawn from a 
long-lived, presumably stellar population, while the group of less active, low-velocity stars represent 
a younger (< 1 Gyr) population . 

Despite the smaller sample size, the properties of L dwarfs are of considerable interest. The 
Sample BC L0-L4 dwarf velocities are typical of an old disk population. Even excluding 2M1300 
due to unusual color, we find atan = 56 km/s, while including it we find atan = 70 km/s. Only two 
of the Sample BC L dwarfs show emission - the low velocity 2M1108 has the strongest emission at 
7.8A, while the high velocity 2M1506 has weak lA emission. While the velocities of only two L 
dwarfs are not definitive, the velocity distribution of the inactive L dwarfs suggest they are mostly 
old. Adding the information provided by the work of K99 and KOO to our data provides strong clues, 
that just as in the ultracool M dwarfs, the traditional activity-age relationship is broken, perhaps 
even reversed. L dwarfs that show lithium absorption are necessarily younger and lower mass than 
L dwarfs of the same spectral type which have destroyed lithium. Thus, using the traditional stellar 
age-activity relation, one would expect them to be more chromospherically active. Consider the LI 
to L4.5 dwarfs, where lithium is detectable even at the low resolution of the K99/K00 Keck LRIS 
observations. Only one L dwarf, Kelu-1, shows both Ha emission and lithium absorption.'^ Eleven 
other such L dwarfs show Ha emission but do not have lithium absorption. Twelve L dwarfs show 
lithium absorption but do not have Ha emission (four of these have marginal Ha detections or noise 
consistent with emission of less than 2A). While many L1-L4 dwarfs have neither Ha emission nor 
lithium absorption, it seems clear that the chromospherically active L dwarfs are drawn from an 
older, more massive population than the lithium L dwarfs. Beyond L4.5, there are no definite cases 
of Ha emission, although lithium absorption is present for ~ 50% of the L dwarfs. 

Brown dwarfs have also been identified in nearby young clusters. In Figure 9, we compare the 



'^It is interesting to note that Basri (1999) finds that Kelu-1 is rotating extremely rapidly: 80 km/s. 
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activity of our field dwarfs to young brown dwarfs. Shown are young brown dwarfs from the a Ori 
cluster (Bejar, Zapatcro Osorio & Rebolo 1999; Zapatero Osorio et al. 1999) and the young brown 
in the p Oph cloud (Luhman, Liebert &: Rieke 1997) - both these clusters arc less than 10 Myr 
old. Also shown are confirmed Plieades brown dwarfs (Martin, Rebolo, & Zapatero-Osorio 1996; 
Zapatero Osorio et al. 1997; Martin et al. 1998a,b) with age ~ 10^ years. In order to suggest the 
age evolution of the field M dwarfs, those with vtan < 20 km/s have been marked as open symbols. 
These low-velocity dwarfs arc likely to be younger than the other field dwarfs. It is evident that the 
Plieades M8 and later brown dwarfs are not more active that the typical field dwarfs, although the 
M6-M7 brown dwarfs appear to be more active. Like the young field L dwarfs that have lithium 
aborption, the Pleiades and a Ori L dwarfs do not show emission, even though some older field L 
dwarfs do. In the case of the p Oph brown dwarf, Luhman, Liebert & Rieke (1997) have shown that 
the emission is probably due to accretion from the circumstellar disk or envelope detected in the 
mid-IR. Similar accretion may account for the a Ori strong emitters, while the absence of accretion 
would account for the weak emission in the other half of the a Ori brown dwarf sample. If the 
emission is chromosphcric, then some other factor (such as rotation) is needed to explain the large 
spread in activity levels. It should be noted that these young brown dwarfs are probably hotter 
at a given spectral type (Luhman, Liebert & Rieke 1997), which suggests that if temperature were 
plotted the brown dwarfs would appear even less active compared to field dwarfs. 

We thus summarize the observations: 

1. Although the fraction of dwarfs showing Ha emission reaches 100% at spectral type M7, the 
fraction that show chromospheric activity drops rapidly for later spectral types. 

2. The fraction of energy in chromospheric Ha for those dwarfs that are active drops rapidly as 
a function of spectral type beyond M6. 

3. Low velocity, kinematically young M8.0-M9.5 dwarfs have weaker activity than many higher 
velocity, old M8.0-M9.5 dwarfs 

4. Flaring is common among the M7-M9.5 dwarfs 

5. The early L (L0-L4) dwarfs in Sample BC have old kinematics. 

6. Early L (L1-L4.5) dwarfs with Ha emission are old and massive enough to have burned lithium 

7. L1-L4.5 dwarfs with lithium are unlikely to have Ha emission. 

8. None of the L dwarfs later than L4.5 have Ha emission but half have lithium. 

9. The two known L dwarfs in young clusters do not show Ha emission 

10. Young Pleiades M8-M9 dwarfs are less active than the higher velocity field M dwarfs. 
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6.3. Discussion 

How can these observations be explained? We believe they imply that the maximum activity 
level is a strong function of temperature beyond spectral type M7, with lower temperature objects 
able to maintain less emission. Additionally, beyond spectral type M7 substellar dwarfs tend to 
have less activity than stellar dwarfs. 

In Figure 10, we plot evolutionary sequences from Burrows et al. (1993) and Baraffe et al. 
(1998). In the Baraffe et al. (1998) models the hydrogen burning limit is at 0.072 Mq and the 
lithium burning limit is at O.OSSMq. Also shown is an estimated temperature scale from Reid 
et al. (1999a) based on the arguments made in K99. The models indicate that it takes > 10^ 
years for stars near the hydrogen burning limit to settle into the MB and cooler temperatures 
(Figure 10). M8 and cooler temperatures are possible at younger ages, but they are substellar 
objects which continue to cool with time. Thus, the low velocity, low activity population is likely 
to be a population of substellar /transition objects. By the time they are older than 1 Gyr, they 
appear as L dwarfs or even cooler T dwarfs. Thus, the comparison of low velocity M8-M9.5 dwarfs 
to high velocity M8-M9.5 dwarfs is the same as a comparison of younger, lower mass O.OTMq) 
objects to older, higher mass objects (~ O.OSMq) at the same temperature. The same occurs when 
comparing the K99/K00 lithium L1-L4 (~ O.OSSM©) dwarfs to the K99/K00 non-lithium L1-L4 
(~ O.O7M0) dwarfs. These mass estimates are only meant as illustrative values - mass estimates 
are subject to many uncertainties and a range of masses and ages will be sampled. 

Both the field M8-M9 dwarf and field L dwarf observations show that at a given spectral type, 
the less massive dwarfs are less active, even though they are younger. It is perhaps worth noting 
that theoretical models suggest that the lower mass objects will be slightly more luminous with a 
smaller surface gravity (Burrows et al. 1993; Baraffe ct al. 1998). The comparison to cluster brown 
dwarfs is consistent with this trend (Figure 9). The Plieades brown dwarfs which have cooled to 
these ultracool M and L temperatures are less active than old field ultracool M and L dwarfs — but 
more active than the low velocity ultracool M dwarfs. There are thus strong suggestions, as already 
noted in Section 6.1, that activity levels do decrease with age in brown dwarfs, and that therefore 
activity levels are dependent upon temperature, mass, and age. The importance of accretion needs 
to be investigated for the youngest (< 10 Myr) ages. The dominant effect is temperature, as both 
young and old objects show the rapid fall in activity levels beyond spectral type M7. 

At the same time, the models suggest that stars, or at least very-long lived hydrogen burning 
transition objects, are likely to exist down to LO - L4 temperatures. This is completely consistent 
with our empirical finding that the early L dwarfs have old kinematics. We note that Kirkpatrick 
et al. (1999a) find a temperature of 1900 it 100 K for the L4 dwarf GD 165B and constrain the age 
to be greater than 1.2 Gyr using updated models and the white dwarf primary's cooling age and 
argue it is just below the substellar limit, near the transition region between stars and brown dwarfs 
(formally, they actually derive the minimum stellar mass using the models). L dwarfs with lithium 
must be below the lithium burning limit (~ O.OSSM©; Chabrier &: Baraffe 1997) and younger than 
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1 Gyr (Figure 10). The inactivity of these lithium L dwarfs demonstrates that the lowest mass 
objects cannot sustain significant activity at a temperature (or luminosity) that is adequate for 
sustaining some activity in older but more massive dwarfs. 

The decline in the frequency of activity may be associated with two effects. First, as later 
spectral types are considered, a larger fraction of very low mass (substcUar) objects contribute, and 
these are more likely to be inactive in field samples. Second, activity among the L dwarfs may die 
out in time, since the observed high velocity L dwarfs are inactive - although we cannot tell if they 
were ever active. It would be of great interest to find whether or not the early K99/K00 L dwarfs 
which are active have low or high velocities. The high-velocity, low-activity M9.5 dwarf 2MASSW 
J0109216-I-294925 may be a young brown dwarf that happens to have high velocity, an old object 
which has never been active, or an old stellar M9 dwarf whose chromospheric activity has declined 
with age. In any case, it is worthy of additional study. We note that the high observed flare rate 
implies that the rotation rate may decrease with time, even among the "inactive, rapid rotators" 
if the flaring is associated with mass loss and/or a stellar wind. We speculate this may provide a 
mechanism for the evolution of activity. Additional observations are needed to determine what the 
rotational velocities are as a function of mass, spectral type, and age. 

What fraction of the ultracool M dwarfs are likely to be substellar? While we cannot identify 
which individual objects are brown dwarfs, we can identify a number of probably young objects. 
Three of our M8.0-M9.5 dwarfs show no Ha emission and very low velocity; another three have 
equally low velocities and EWna < 3A. Out of a total population of 32 M8.0-M9.5 dwarfs, our 
data suggest that 10 — 20% are brown dwarfs. These objects should be more likely to have lithium 
absorption (like LP 944-20) , but most of the brown dwarfs will be massive enough to burn lithium. 
Indeed, we note that none show the distinctive signs of low surface gravity that characterize the 
Pleiades M8 brown dwarfs Teide 1 and Calar 3 (Martin, Rebolo, & Zapatero-Osorio 1996), so none 
of our targets arc very young {r^ 10^ years). Approximately ten objects belong to the low velocity, 
low activity group in the lower left of Figure 8 - that is a third of the sample, but some fraction of 
these will stabilize as hydrogen burning L dwarfs, in order to account for the observed popoulation 
of high velocity early L dwarfs. These fractions will be somewhat overestimated for the Galactic 
disk population, since the old, large scale height population will be imdcrrepresented locally. Other 
effects may also be important, such as the fact that wc have estimated distances for all dwarfs 
using one color-absolute magnitude relation. Adding kinematic ages, as in this study, provides an 
additional constraint on the modelling necessary to determine the field substellar mass function 
(Reid et al. 1999a). 

The nature of the ultracool M dwarfs has been debated for some time in the literature. While 
convential wisdom suggested that most if not all field ultracool M dwarfs are stellar, many sug- 
gestions that most ultracool M dwarfs are substellar have been made, most of which have been 
discredited. We remark that Bcssell (1991) noted that the ultracool M dwarfs are expected to be 
a mixture of young brown dwarfs and older stars - and he also noted that there was a paucity of 
high proper motion ultracool M dwarfs expected from the stellar population in the LHS catalogue 
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(Luyten 1979). Our study has identified a number of high proper M dwarfs which appear on the red 
POSS plates but were overlooked for the LHS catalogue - evidently, the faintest of these targets on 
the blue plate precluded their detection by Luytcn and contributed to the effect noted by Bessell. 
Our results show that most ultracool M dwarfs are old, and hence stellar, but perhaps 10-20% are 
a younger population of brown dwarfs. 

We end our discussion with a few caveats due to our photometric selection. Both the relative 
numbers and kinematics of "active" and "inactive" M dwarfs will be changed if one group is 
preferentially brighter at Mk for its J — K color. Indeed, the inactive brown dwarf LP 944-20 lies 
one magnitude below the active M dwarf LHS 2397a in the J — K, Mk HR diagram. Preliminary 
USNO parallaxes show instrinsic dispersion in the HR diagram for the late M and L dwarfs (Dahn, 
private communication). If inactive ultracool M dwarfs are sub luminous compared to our adopted 
relation, we will have overestimated their velocities; correspondingly, if the more active dwarfs are 
"superluminous," we have underestimated their velocities. Fortunately, this would only strengthen 
our evidence that active ultracool M dwarfs are older. The intrinsic dispersion presumably depends 
upon such ill-understood factors as metallicity, age, surface gravity, and dust formation. Another 
possible bias on activity levels is that we favor the inclusion of unresolved binaries. Amongst the 
earlier M dwarfs, very short period systems have enhanced chromospheric activity due to tidal 
effects maintaining high rotation rates(Young, Sadjadi, &: Harlan 1987) — however, even if this 
mechanism works in the ultracool M dwarfs, which seems unlikely if the rotation-activity relation 
has broken down, only ~ 5% of earlier type M dwarfs show emission due to this effect, so it should 
be negligible. 



7. Summary 

We show that a sample of bright, nearby ultracool M and L dwarfs can be selected without 
proper motion bias using 2MASS and PMM scans. Our initial samples include high proper motion 
objects, visible on the POSS plates, that should be added to an updated version of the LHS 
Catalogue, and one M8.0 dwarf with a photometric parallax that places it within 10 parsecs. We 
intend to continue this study in order to complete the nearby star catalog for the lowest mass stars. 

Using our initial sample, we estimate the space density of dwarfs near the hydrogen-burning 
limit. We show that the dropoff near the hydrogen burning limit in the five and eight parsec nearby 
star samples is likely to be due to incompleteness. This is more consistent with a smooth relation 
across the hydrogen burning limit. Trigonometric parallaxes and searches for companions will help 
improve the space density estimate. 

Most importantly, we use our spectroscopic observations of our well-defined sample to explore 
the relationships between age, kinematics, and chromospheric activity for the ultracool M and 
L dwarfs. We show that the observations can be understood if activity is primarily related to 
temperature and secondarily mass and age, and that lower mass (substellar) objects have weaker 
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chromosperes. Thus, the classical relation that strong Ha emission implies youth is not valid for 
these dwarfs. Instead, strong Ha emitters in the field are likely to be old (> 1 Gyr) stars, while 
weaker emitters are often young (< 1 Gyr), lower-mass brown dwarfs. This does not exclude the 
idea that for a given dwarf. Ha activity declines with age - but spectral type (temperature) is the 
observable in the field. The local population of ultracool M dwarfs apparently consists both of the 
most massive (lithium burning) brown dwarfs and the lowest mass (hydrogen burning) stars, with 
the substellar objects making up a significant fraction of the sample. The early L (L0-L4) dwarfs 
are consistent with an old, at least partially stellar population. The evidence thus suggests, as do 
some models, that early L dwarfs can be stable hydrogen-burning stars. Expansion of the sample 
with follow-up observations should clarify the relative contribution of stars and brown dwarfs to 
these temperature ranges. 
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,5 


4, 


.9 


C 


2MASSW J2313472+211729"' 


23:13:47.29 


+21:17:29, 


,5 


11, 


.43 


10 


.75 


10 


.42 


M6 




6, 


.0 


B 


2MASSW J2331016-040618 


23:31:01.63 


-04:06:18. 


,6 


12, 


.94 


12, 


.29 


11 


.93 


M8. 


,0 


5, 


.4 


C 


2MASSI J2334394+ 193304 


23:34:39.44 


+ 19:33:04. 


,2 


12, 


.77 


12, 


.07 


11, 


.64 


M8. 


,0 


22, 


.6 


B 


2MASSI J2336439+215338'^ 


23:36:43.92 


+21:53:38. 


,7 


12, 


.76 


12, 


.10 


11 


.71 


M7. 


,0 


7, 


.7 


B 


2MASSW J2347368+270206 


23:47:36.80 


+27:02:06. 


,8 


13, 


.19 


12, 


.45 


12, 


.00 


M9. 


,0 


3 


.0 


B 


2MASSW J2349489+122438° 


23:49:48.99 


+12:24:38. 


,8 


12, 


.62 


11, 


.95 


11, 


.56 


M8. 


,0 


3 


.5 


C 


2MASSW J2358290+270205P 


23:58:29.00 


+27:02:05. 


,5 


12, 


.71 


12, 


.05 


11 


.68 


M6 




5, 


.9 


B 



^LP 349-25 
''LP 356-770 
'^LP 413-53 
^LP 427-38 
^LP 326-21 
^LP 328-36 
SLP 44-162 
'^LP 872-22 
^LP 759-17 
JLP 288-31 
•^LP 460-44 
'LP 345-18 
"^LP 461-11 
°LP 402-58 
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Table 2. L Dwarf Data 



Name 


CrH-a 


Rb-b/TiO-b 


Cs-a/VO-b 


K I fit 


Type 


2MASSW J1108307+683017 


1.27(0-1) 


0.88(1) 


0.84(1) 


(0) 


LI V 


2MASSW J1300425+191235 


1.53(2) 


0.81(1) 


0.81(0-1) 


(1) 


LI V 


2MASSW J1506544+132106 


1.44(1-2) 


1.18(3) 


1.13(3) 


(3) 


L3 V 


2MASSW J1658037+ 702701 


1.26(0-1) 


0.79(1) 


0.81(0-1) 


(2) 


LI V 
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Table 3. Kinematics and Activity 



Name 




fJ-a 


1^5 




log 4^ 


BR 1222-1221 




.\J 


-0 


9fi9 


-0 


1 87 


25 


-4.70 


LHS 2243 


1 ft 

±\J. 


ft 


n 
-u. 







.^U J- 


38 


-3.57 


LHS 2397a 


19 







50Q 


-0 


081 


29 


-4.22 


LP 412-31 


1 1 

± A. . 


7 


n 


34Q 


-0 


951 


24 


-4.45 


LP 415-20 


22, 




n 


127 


-0 


036 


14 


-4.56 


LP 475-855 


99 


9 


n 




n 

u. 


.uxu 


11 


-3.23 


2MASSW J0027559-F221932 


8 

o. 


. (J 


n 


40 


-0 


1 79 


17 


-4.52 


2MASSW J0036159-hl82110 


1 1 

± J. , 


1 

. J. 


n 
u. 


.oo ( 




1 04 


44 




2MASSW J0109216-F294925 


18 

-LO. 


7 


1 







348 


95 




2MASSW J0140026+270150 


1 Q 


A 




u. 


Ofil 


n 

-u. 


91^9 


24 




2MASSI J0149089+295613 


1 7 

J. 1 . 


4 





907 


-0 


466 


42 


-4.62 


2MASSI J0240295+283257 


99 


7 


n 


04fi 


-0 


1 99 


21 


-4.27 


2MASSI J0253202-F271333 


90 


1 

. ± 


n 


■^70 


n 


088 


36 


-4.31 


2MASSI J0330050-F240528 


90 


1 

. J. 


n 


18^ 


-0 


03Q 


18 


-3.91 


2MASSI J0335020-F234235 


IQ 


9 


n 


058 


-0 


043 


7 


-4.71 


2MASSW J0350573+181806 


IQ 


q 


n 

u. 


18Q 


-0 

u. 


04Q 


18 




2MASSW J0354013+231633 


99 


.o 


-0 


lfi8 


n 


064 


19 


-4.67 


2MASSI J0746425-h200032 


"1 


.4: 


n 


■ OOO 


n 

u. 




18 




2MASSI J0810586-M42039 


90 


Q 
.O 


n 




n 

u. 


1 98 


13 


-4.69 


2MASSI J0818580-F233352 


1 7 


Q 


n 


97^ 


n 

u. 




35 


-4.31 


2MASSI J0925348-M 70441 


ZZ. 


r 
.0 


-U. 


.Z6Z 


U, 


.UiU 


25 


-4.28 


2MASSW J0952219-192431 


15, 


.4 


-0. 


,077 


-0, 


.104 


9 


-3.96 


2MASSI J1024099+181553 


18, 


.2 


-0. 


,144 


-0, 


.070 


14 


-4.73 


2MASSW J1108307-F683017 


12, 


.8 


-0. 


,226 


-0, 


.194 


18 


-5.44 


2MASSW J1127534-F741107 


24, 


.3 


-0. 


,016 


-0, 


.030 


4 


-4.84 


2MASSW J1200329-F204851 


24, 


.1 


-0. 


,159 


0, 


.232 


32 


-4.44 


2MASSW J1300425+191235 


20, 


.3 


-0. 


,789 


-1, 


.238 


141 




2MASSW J1311391-F803222 


21, 


.3 


-0. 


,068 


-0, 


.348 


36 


-4.86 


2MASSW J1336504+475131 


22, 


.5 


0, 


,111 


-0, 


.016 


12 


-4.47 


2MASSW J1344582+771551 


23, 


.7 


0. 


,072 


-0, 


.005 


8 


-4.68 


2MASSW J1403223-F300755 


21, 


.4 


-0, 


,788 


0, 


.042 


80 


-4.07 


2MASSW J1421314-M82740 


19, 


.6 


-0. 


,744 


-0, 


.182 


71 


-4.97 


2MASSI J1426316-M55701 


20, 


.0 


0. 


,108 


-0, 


.056 


12 


-5.40 
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Table 3 — Continued 



Name 


dphot 


fJ-a 


IJ-S 




log^ 


2MASSW J1439283+192915 


18, 


.5 


-1. 


,245 


0, 


.392 


114 




2MASSW J1444171+300214 


12, 


.5 


-0. 


,101 


-0, 


.336 


21 


-4.42 


2MASSW J1457396+451716 


20, 


.7 


-0. 


,191 


0, 


.100 


21 


-4.76 


2MASSW J1506544+132106 


12, 


.1 


-1. 


,092 


0, 


.001 


63 


-5.72 


2MASSW J 1546054+374946 


20, 


.2 


-0. 


,020 


-0, 


.120 


12 


-3.90 


2MASSW J1550381+304103 


24, 


.2 


-0, 


,112 


0, 


,107 


18 


-3.73 


2MASSW J1551066+645704 


20, 


,6 


-0, 


,220 


0, 


.010 


22 


-4.39 


2MASSW J1553199+140033 


21, 


.1 


-0. 


,659 


0, 


.072 


66 


-4.48 


2MASSW J1627279+810507 


21, 


.3 


-0. 


,209 


0, 


.338 


40 


-4.77 


2MASSW J1635192+422305 


22, 


.4 


-0. 


,073 


-0, 


.010 


8 


-5.12 


2MASSW J1658037+702701 


17, 


.4 


-0. 


,136 


-0, 


.315 


28 




2MASSW J1707183+643933 


17, 


,1 


0, 


,226 


-0, 


.091 


20 


-4.49 


2MASSW J1733189+463359 


17, 


.6 


0, 


,044 


-0, 


.257 


22 


-5.27 


2MASSW J1750129+442404 


23, 


.4 


-0, 


,018 


0, 


.151 


17 


-4.87 


2MASSW J1757154+704201 


11, 


.7 


0. 


,006 


0, 


.338 


19 


-4.83 


2MASSW J2049197-194432 


21, 


.9 


0. 


,193 


-0, 


.260 


34 


-4.24 


2MASSW J2140293+162518 


20, 


.7 


-0, 


,008 


-0, 


.102 


10 




2MASSW J2 147446-264406 


22, 


.9 


-0, 


,054 


-0, 


.232 


26 


-5.04 


2MASSW J2206228-204705 


18, 


,4 


0, 


,001 


-0, 


.065 


6 


-4.59 


2MASSI J2234138+235956 


17, 


.6 


0, 


,829 


-0, 


.034 


69 


-4.91 


2MASSI J2235490+184029 


17, 


.3 


0. 


,326 


0, 


.042 


27 


-4.41 


2MASSW J2306292-050227 


11, 


.3 


0. 


,889 


-0, 


.420 


53 


-4.61 


2MASSW J2331016-040618 


26, 


.3 


0, 


,401 


-0, 


.231 


58 


-4.72 


2MASSI J2334394+ 193304 


20, 


.0 


-0, 


,236 


-0, 


.117 


25 


-4.10 


2MASSI J2336439+215338 


22, 


.4 


0, 


,379 


0, 


.024 


40 


-4.29 


2MASSW J2347368+270206 


22, 


.2 


0, 


,313 


0, 


.033 


33 


-5.21 


2MASSW J2349489+122438 


20, 


.7 


0. 


,025 


-0, 


.189 


19 


-4.78 
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Table 4. Space Densities 



Sp. Type 


N 


$ 








units 




</ > 

Vmax 




M8.0-M8.5 


17. 


1.90 


0.47 


IQ- 


-3 


stars pc" 


-3 


0.56 


0.07 


M9.0-M9.5 


15. 


2.57 


0.69 


10- 


-3 


stars pc- 


-3 


0.71 


0.07 


M8.0-M9.5 


32. 


4.46 


0.83 


10" 


-3 


stars pc- 


-3 


0.63 


0.05 


L0.0-L4.5 


7. 


2.11 


0.92 


10" 


-3 


stars pc" 


-3 


0.53 


0.11 


M8.0-L4.5 


39. 


6.57 


1.24 


10" 


-3 


stars pc" 


-3 


0.61 


0.05 


M8.0-M9.5 


32. 


4.75 


0.89 


10" 


-3 


stars pc" 


-3 —1 

mag 


0.63 


0.05 


M8.0-L4.5 


39. 


4.38 


0.83 


10- 


-3 


stars pc" 


—1 

mag 


0.61 


0.05 



Table 5. Flares 



Name 


Sp. 


Ha 


Source 


Flare Ha 


Source 


2MASSW J0149089+295613 


M9.5 


11.0 


L99 


300 


L99 


2MASSW J2234138+235956 


M9.5 


4.4 


KPNO 


20 


Keck 


LHS 2397a 


M8 


15 


LCO 


? 


B91 


LHS 2243 


M8 


1.3 


M94 


44 


LCO 


LP 475-855 


M7 


7 


Keck 


40.5 


LCO 



Note. — Sources are: LCO (Las Campanas, this paper), KPNO (Kitt 
Peak, this paper), Keck (Keck, this paper), B91 (Bessell 1991), M94 (Martin, 
Rebolo, & Maguzzu 1994), L99 (Liebert et al. 1999) 
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Fig. 1. — Absolute magnitudes for M7 and later dwarfs (solid points) with the linear fit Mk = 
7.593 + 2.25 x ( J - Kg) shown. The scatter about this fit is cr = 0.36 magnitudes. 2MASS data 
for Hyades members from Gizis et al. (1999) is also shown to illustrate the steepness of of main 
sequence for for M0-M7 dwarfs, which may lead to errors in the distance estimates for M7 dwarfs. 



-31 - 




Fig. 2. — Color-color diagram for the December sample using the photographic R and 2MASS J 
and Kg magnitudes. M8 and later dwarfs are solid circles, M7 dwarfs are open squares, and M6 
and earlier dwarfs are crosses. A selection on R-K^ > 4.9 would improve our efficiency without 
excluding the ultracool M dwarfs. 
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Fig. 3. — 2MASS J-K^ color as a function of spectral type. Each bin is 0.05 magnitudes wide, 
approximately equal to the 2MASS uncertainty. Note the correlation between infrared color and 
far-red spectral type. 
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Fig. 4. — 2MASS near-infrared color-color diagram. The solid lines indicate our selection criterion. 
Note that all the M and L dwarfs lie well below the line intended to exlude giants. 
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Fig. 5. — Our observed space density of cool (M8.0-L4.5) dwarfs compared to the the Reid & Gizis 
(1997) eight parsec sample (as updated in Reid et al. 1999). The solid histogram count known 
secondaries, while the dotted histogram excludes them. The steep dropoff at Mk > 10.0 seen in 
both the eight and 5.2 parsec samples is moderated by our space density. 
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Fig. 6. — The observed percentage of Ha emission line dwarfs amongst K5 to M6.5 dwarfs (HGR, 
open squares), M6 to M9.5 dwarfs (this paper, soHd circles), and LO to L8 dwarfs (K99;K00, 
crosses). K5 and K7 dwarfs are plotted as -2 and -1 respectively, while the L dwarfs are plotted 
with 10 added to their subclass. The HGR M6 and M6.5 dwarfs may be affected by kinematic 
bias, leading to an underestimate of the number of emission stars. The solid line connects the three 
studies, adopting this paper's values for M6 dwarfs over the HGR values. 
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Fig. 7. — The Ha luminosity relative to the bolometric luminosity as a function of spectral type 
for our ultracool M dwarfs (solid circles), the earlier M dwarfs of HGR (open squares), and the L 
dwarfs of K99 (crosses). For our ultracool M dwarfs, approximate upper limits are plotted assuming 
an Ha equivalent width of 2 A. The dotted line at -3.9 is the level at which any M dwarf would be 
observed in emission. None of the M8 or later dwarfs have activity levels above the -3.9 level. 
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Fig. 8. — Ha equivalent width as a function of tangential velocity. M8.0 - MS. 5 dwarfs are plotted 
as stars and M9.0 - M9.5 dwarfs are plotted as solid circles. All M8-M9 dwarfs with Ha equiva- 
lent widths over lOAhave large space velocities. Note 2MASSW J0109216+294925, which has no 
emission but a high tangential velocity. 



-38- 



• • • CO 



lOO o 



I 



□ * 



o 
e o 



o 



o o 



•oo o 

• • •••• 1 

• • 9 

• • □•••c* * 

• (2) • O 

• O ©O" o • 

• : • |» 



Fig. 9. — The Ha luminosity relative to the bolometric luminosity as a function of spectral type 
for both cluster brown dwarfs and field dwarfs. Brown dwarfs from the a Ori cluster (< 10'' years), 
p Oph (< 10^ years), Pleiades (~ 10^ years) are shown as open squares, open triangles, and six- 
pointed stars respectively. Note that both cluster L dwarfs have only upper limits on the detected 
Ha emission. The field M dwarfs are plotted as open circles if vtan < 20 km/s and solid circles for 
higher velocities. 
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Fig. 10. — Model calculations of brown dwarfs and the lowest mass stars by Baraffe et al. (1998), 
including preliminary models for lower masses and younger ages (Baraffe & Chabrier, priv. comm.) 
and for Burrows et al. (1993, 1997). Along the left axis, the estimated temperature scale of K99 
and Raid et al. (1999) is indicated. The hydrogen burning limit is O.O72M0 and the lithium burning 
limit is 0.055Mq (Chabrier k. Baraffe 1997). The models suggest that our spectral range will be 
populated by stars with age > 1 Gyr, transition brown dwarfs burning lithium with 0.4 < age < 1 
Gyr, and brown dwarfs with lithium with age < 0.4 Gyr. Both the model temperatures and 
the spectral-type temperatures are uncertain. These Baraffe et al. models use grainless model 
atmospheres. 



